clear all
clc
%%%%% load data

MRIO2016=xlsread('Russia MRIO 2016');

A=MRIO2016(1:5015,1:5015)./Output';
F=MRIO2016(1:5015,5016:(5015+7*85));
E=MRIO2016(5020,1:5015);

E(isnan(E))=0;
E(isinf(E))=0;
A(isnan(A))=0;
A(isinf(A))=0;

I=eye(size(A));
L=(I-A)^-1;

%%%%% calculate

% 
for i=1:85
    F1(:,i)=sum(F(:,(i-1)*7+1:(i-1)*7+7),2);
    CF2{i}=diag(E)*L*diag(F1(:,i));
end


% 
for i=1:85
    CF3=CF2{i}; 
    CF4= sum(CF3,2); 
    CF4_1(:,i)= sum(CF3,2);
    for j=1:85
      CF5(j,i)=sum(CF4((j-1)*59+1:(j-1)*59+59,:),1);
    end
end


% 
for i=1:85
    CF3=CF2{i};
    CF41= sum(CF3,1);
    CF42(i,:)=CF41;
    
    for j=1:85
      CF51(j,i)=sum(CF41(:,(j-1)*59+1:(j-1)*59+59),2);
    end
end
